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Abstract 

We introduce a novel approach, a Dense Shell Method (DSM), for measuring distances for cosmology. 
It is based on original Baade idea to relate absolute difference of photospheric radii with photospheric 
velocity. We demonstrate that this idea works: the new method does not rely on the Cosmic Distance 
Ladder and gives satisfactory results for the most luminous Type Iln Supernovae. This allows one to 
make them good primary distance indicators for cosmology. Fixing correction factors for illustration, we 
obtain with this method the median distance of sa 68l}g(68%CL) Mpc to SN 2006gy and median Hubble 
parameter 79±??(68%CL) km/s/Mpc. 

1 Introduction 

Supernovae are among the most luminous phenomena in the Universe, and they can serve as cosmological 
distance indicators. In some cases one can use a standard candle method. Nobel prize 2011 in physics 
is given "for the discovery of the accelerating expansion of the Universe through observations of distant 
supernovae" . Actually, Type la supernovae have been used for this. 

Although SNe la are not uniform in luminosity, they can be standardised. The standardisation is based 
on statistical correlations found for nearby events [T] [5] . Thus they are secondary distance indicators, see 
reviews, e.g. [3j[4]. 

Type II supernovae, on the other hand, have a much larger variance in luminosity and therefore cannot 
provide an accurate distance by photometry alone. Nevertheless, their great advantage is the possibility 
of direct measurement of distance, e.g. by Expanding Photosphere Method (EPM) [5] when applied to 
SNe IIP. The development of EPM is the spectral-fitting expanding atmosphere method (SEAM) 6 . Thus, 
Type II supernovae are interesting because there are ways to make them primary distance indicators. A 
standard candle assumption and its calibration is not needed for direct methods. Applications of SNe IIP 
in cosmography do not depend upon the steps of Cosmic Distance Ladder avoiding their systematic and 
statistical errors. 

Due to absolute weakness of SNe IIP they cannot be used at large cosmological distances. In this paper 
we introduce a novel approach to measuring distances for cosmology with the help of the most luminous 
Type Iln Supernovae. The method is based on the formation of an expanding dense shell in SN Iln and 
allows one to find a linear size of a supernova shell in absolute units and distance to it. This Dense Shell 
Method (DSM) is partly based on ideas introduced in EPM and SEAM, and partly in Expanding Shock 
Front Method (ESM) [7J used for SNR 1993 J. 
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2 Classical Baade-Wesselink and Kirshner-Kwan methods 



All researches using EPM for supernovae cite papers [S] and [5]. Actually, EPM introduced by [5] differs 
from the classical Baade-Wesselink (BW) method. 

Here we repeat briefly the steps of BW approach which we apply in our new method. 

Measuring colour and flux at two different times, t\ and t2, one finds the ratio of the star's radii, R2/R1 
(the same can be found from interferometry). Using weak lines which are believed to be formed near the 
photosphere one can measure, in principle, the photospheric speed « p h = dR p ^/dt. Then J 2 v p ^dt would 
give Aiiph = R2 — Ri- Knowing R2/R1 and R2 — Ri, it is easy to solve for the radii. The ratio of fluxes 
gives the distance D: 

D = Rph \l F^owttdT ■ (1) 

Actually, finding the distance by equation ([T]) with F v (model) is equivalent to Spectral-fitting Expanding 
Atmosphere Method (SEAM) 6 . The original BW method is based on a simplifying assumption of a 
diluted supernova blackbody spectrum, 

F„(model) = -k&B„{T c ) . (2) 

Here the relation of a true photospheric intensity with blackbody brightness B V {T C ) is accounted for by a 
correction factor f„. This factor is often called the dilution factor (a ratio of a thermalisation radius to 
i? p h). Thus, the distance is: 

° = CvR ^\l FjotservSy- (3) 

Apart from the correction for the dilution one needs also a correction for limb darkening, or brightening, 
for the ratio of pulsation velocity to the radial velocity accounted for by projection factor p. See, e.g. 
[1UI 111! 1121 1131 114] and references therein for a discussion of those non-trivial questions on the projection 
factor and other problems related to BW method in Cepheids. 

In reality, one can measure directly only the matter velocity v m on the photospheric level. The as- 
sumption v m = w p h does not work (as a rule) in exploding stars. Even for Cepheids this was questioned 
already by [15J . Velocity of matter at the photosphere of a supernova is not at all dRph/dt. The v p h and 
Vm may even have different signs. That is why the main idea of EPM for SNe is different from BW. 

Kirshner and Kwan \E\ also used the weak lines to measure the matter velocity on photospheric level, 
Vm, but they never put v m equal to dR p h/dt. That is why the EPM for supernovae should be called not 
the Baade-Wesselink method, but more properly the Kirshner-Kwan method (KK). They determine the 
photospheric radius from the relation 

R P h = v p h(t — t ) , (4) 

where to is the constant close to the explosion epoch. This relation is based on the assumption of free 
expansion. If R p h is obtained, the distance D to the supernova is found from equation ©. 



3 Direct distance determination by the new method 

Let us introduce briefly the essence of the new Dense Shell Method (DSM). 

Supernovae of type Iln, contrary to SNe IIP, do not enter the coasting free expansion phase and both 
EPM [5] and SEAM [6] are not directly applicable. Nevertheless, in SN Iln case we can use slightly 
modified classical BW method. There is a lot of dense matter around the supernova and the shock cannot 
break out through the circumstellar shell for months or even years. Yet, it is clear from our results on 
SNe Iln |161 117] that all matter behind the shock is cooled down by radiation and compressed into a cold 
dense shell. One has to measure wide emission components of lines and determine velocity of matter in the 
dense shell v^s (with highest possible accuracy). Since forward and reverse shocks are both glued together 
in this shell the photosphere moves with the matter as well. In the dense shell, v p h is exactly equal to the 
rate of change of R p h, i.e. w p h = dR p h/dt = v^s - and this can be measured. Everything looks as Baade 
suggested already in 1920s! 

First, we formulate the DSM for broad-band flux F and integrated correction factor £. The observed 
flux is F = C, 2 R ph ivB(T)/ D 2 , where B(T) is the blackbody intensity and D is the photometric distance. 
Then yF = (,R p h^J nB(T) / D. The effective blackbody temperature T is measurable, as well as dR p ^ and 
dy/F, while D does not change. 

Hence, if T and ( are almost constant between the two measurements, we have 

dVF = {dR ph ^nB(T)/D, (5) 
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and 

D = C,dR pil ^-KB(T)/d^fF. (6) 

Thus, measuring dvF, dR p ^ and T, and calculating £ from a model, we find the distance D by a direct 

method without any ladder of cosmological distances. 

One may limit oneself with this "two-point" method for quick evaluation of the distance 

However, this quick estimate may result in a large error when fluxes are close to each other and dy/F 

is small in denominator of equation (JSJl . For a more accurate treatment one has to develop the new robust 

technique for monochromatic or broad-band fluxes, correction factors, and variable colour temperature 

T{t,u). 

If temperature T changes significantly with time t and frequency v, we have to rely on the evolution of 
R p h, which is controlled by the changes of radii d_R p h(i) taken from observations. Moreover, we have also 
to use a model to calculate a correction factor and theoretical flux F v . 

Assume that the observations are sufficiently frequent to allow us to measure the increments in radius 
dRph = Vphdt for a number of points, where dt is a difference of time of the successive observations. 

Let the initial radius (unknown to us) is Rq, and Ri = Rq + ARi for i — 1, 2, 3, . . ., where ARi is 
already known from the dR integration over time. 

Then 

Ci(Ro + ARi)\B4T c „i) = 1Q 0AA "D 2 F VI (7) 

or, by taking the root, 

Cui{Ro + AIU)y/irB v (T C ui) = 10°' 2A " Dy/FU. (8) 

Here A v is the extinction in stellar magnitudes for the frequency v. A good model gives us a set of the 
Ci/ii Tcui for all observational points. From the measured F vi , ARi we can find Rq and the combination 
a s = 10 0AAs D 2 (where instead of v we use index s labelling one of the broad-band filters) by the least 
squares method. 

To find the distance D we need to know A s from the astronomical observations, or we can try to get 
it from equation ((8} written for different spectral filters. 
Knowing Rq we obtain the set of equations: 

1U D = a s , 

This gives us difference A 3 \ — A S 2, and with the help of, e.g., [18| law one may find A s . 

Actually, we have a set of different trial models with different trial distances. All unsuitable models 
(which do not reproduce the shapes of time-dependence of fluxes and colours with reasonable accuracy 
when scaled to a proper distance) are discarded. "Suitable" means that they reproduce the observed values 
of velocity v, temperature T, and circumstellar envelope density p. Scaling means that they have different 
radii of the dense shell for any given time. Hence they correspond to different distances to SN. 

The high luminosity of type Iln supernovae is explained by inflowing matter merging with the dense 
shell in highly radiative shocks. From the continuity of mass we find: 

'' s " I (9) 



VS — Vi po 

where vs is the shock velocity, po is the density ahead of the front, and vi, pi are velocity and density 
behind the front respectively, see Fig. 16 in [TS] and Figs. S10, Sll in |17| . Due to the extremely high 
density contrast we have vi ~ vs (hence merging of forward and reverse shocks). This pattern of the flow 
is obtained when all kinetic energy of inflowing matter in the rest frame of the shock front is radiated 
away: 

^ = = p^? (1Q) 

Thus the observed broad line components permit to determine vs, and the location of of the radiation 
flux creation (and hence the thermalisation radius). As found in our computations, the latter is close to 
the photosphere. This means that the values of £ do not vary too much from model to model and do 
not strongly depend on the photospheric radius. That is, all our "suitable" models will give us a correct 
distances D from solving the system of equations |[8]| in one step despite of the differences in the initially 
assumed distances. 

Other results are obtained by applying EPM of 5 (KK method) to classical type for SNe IIP. First 
[191 120] found that £ depends primarily on temperature T and can be applied for different photospheric 
radii. But this proved to be wrong! E.g. |21j have corrected this statement in case of SN 1999em, where 
they have obtained systematically larger correction factors than [191 120] , Our models for SN 1999em [22] 
support this conclusion quite independently |23| . The correction factor for SNe IIP is more sensitive to the 
model photospheric radius, than in our SNe Iln models with their dense shell: larger is the radius of the 
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Figure 1: Monte Carlo resampling simulation of the distance D to SN 2006gy by DSM method. The 
isocontours of probability distribution function (pdf ) are shown with equal step in pdf. The observations from 
(23] [55] have been used for six different time points from Tabled] (the first point at t = 36.03 d is discarded). 

SNe IIP higher is the scattering in its atmosphere, and hence larger is R p ^ relative to the thermalisation 
radius. Thus, £ depends not only on T, but also on R p h of supernova of type IIP. In our case, we have a 
formation of a dense shell and £ is virtually independent of its radius. Therefore, in EPM or better to say 
in KK method one has to iterate a set of models with system J3J, Q to obtain self-consistent values of £ 
and distance D. 

4 Distance and Hubble's parameter 

For illustration we have taken observational data for SN 2006gy from |24l I25| . Unfortunately, the number 
of epochs for measured temperature is less then the number of observations for fluxes. We have collected 
the suitable data points in Table [1] with interpolation in temperature. 

We have adopted v — 5200 ± 320 km/s from [25]. That is the value corresponding to the rising part of 
the light curve when £ « 1 and the shell does not fragment. 

There are several suggested values for the extinction Ar 26, 24 . We have taken Ar = 1.3 ± 0.25 mag 
following [2B], see discussion in |27j . 

To estimate the confidence intervals of the distance and -Ho we have done a resampling Monte Carlo 
(MC) simulation based on these data. We resampled the values of T, of the stellar magnitude rtiR in 
standard filter R [55], the reddening Ar, and velocity v each with normal distribution having standard 
deviations a from Table [T] 

For obtaining the confidence intervals for the mean and median it was sufficient to do 10 MC tests. 
The plot in Fig. Q] is built with 10 7 samples to obtain a better statistics near the top of the distribution. 

Using all 7 points in Table [1] we have obtained the mean distance D as 63.5 Mpc, and median D w 
62.6 Mpc with 68% confidence interval (-16, +19) Mpc. 

This simulation used the correction factor f = 1, which is close to the values of £ with accuracy about 
~ 10 % found in our radiation hydro models [17] for the growing part of the light curve. Of course, an 
accurate modelling requires building a hydrodynamical model not only for the light curve but also for 
spectral line profiles with account of dilution and projection effects as is being done for recent Cepheid 
models. 

We see that gt of the first point at t — 36.03 d in Table [T] is very large. If we artificially reduce it by a 
factor of 10 we find the mean distance D w 73.4 Mpc, and median D w 72.3 Mpc with the 68% confidence 
interval (-15, +18) Mpc. 

If we discard this point completely we get the mean distance D « 70.3 Mpc, and median D ~ 68.2 Mpc 
with 68% confidence interval (—15, +19) Mpc. Those experiments show that the results are rather robust 
given the level of accuracy of data and models. 
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Table 1: Observations of SN 2006gy 



time, d T, 10 3 K err.(T) ~m~ R crr.(m fl ) 



36.03 12 3 14.72 .03 

40.95 12 .8 14.62 .03 

47.97 12 .8 14.42 .03 

59.92 12 .8 14.27 .03 

71.0, 11 .7 14.22 .03 

82.92 9 .8 14.28 .03 

94.88 8.8 .4 14.49 .03 



This value of distance D is in good agreement with a generally accepted value 71 Mpc, see Fig. [TJ The 
largest thick-line contour in Fig. [1] is about one standard deviation. The error of our value is quite high 
mainly because of the uncertainty of temperature T and reddening Ar. Nevertheless, even this accuracy is 
enough to make quite implausible the suggestion [29] to put SN 2006gy much closer to us, around 10 Mpc. 
The supernova itself "tells" us, that its distance is an order of magnitude larger than 10 Mpc. 

Using the redshift z = 0.0179 for the galaxy NGC 1260, where SN 2006gy has exploded, we obtain the 
Hubble parameter. We do not use directly the formula 



since, e.g., median(_ffo) is not equal cz/median(Z)). So, our values of Ho are computed as a result of MC 
for each individual D. 

Thus, with all data of Table [T] we get the mean Ho ~ 95.2 km/s/Mpc and the median Ho ~ 
85.7 km/s/Mpc with the 68% confidence interval (-20, +29) km/s/Mpc. 

If we reduce the error of the first point by a factor of 10 we find the mean Ho ~ 76.5 km/s/Mpc and 
the median Ho ~ 74.2 km/s/Mpc with 68% confidence interval (—15, +19) km/s/Mpc. 

When the first point is discarded, we get the mean Ho « 81.5 km/s/Mpc and the median Ho ~ 
78.7 km/s/Mpc with 68% confidence interval (-17, +23) km/s/Mpc. 

The latter result is the most reliable, so the median for the Hubble parameter is 



The accuracy is about 30%, mainly influenced by the error in the temperature and the interstellar 
extinction [25] . Of course, this accuracy of Ho is low compared with the one already achieved by other 
techniques, but our value is obtained by the new direct method and does not rely on the Cosmic Distance 
Ladder. 

Statistics of similar objects with more precise reddening can significantly improve the Ho accuracy in 
the future. It is needed to investigate the role of variations of the correction factors in different SN 2006gy 
models to check the robustness of our results. We present here the values for D and Ho only for the 
illustration of the efficiency of the method. 

5 Conclusions 

Now, we can summarise essential features of the new method, DSM (Dense Shell Method), for finding 
cosmological distances with the help of SNe Iln. The method is based on the following steps: 

• Measurement of wide emission components of lines and determination of the velocity at photosphere 
level v m = v p h (with highest possible accuracy). 

• Measurement of narrow components of spectral lines for estimating properties (density, velocity) of 
circumstellar envelope. One does not need a very high accuracy of measurements and modelling here. 

• Building of a set of best fitting models ("suitable") for broad band photometry and speed v p h, for a 
set of trial distances satisfying the constraints for the circumstellar envelope found from narrow lines. 

• Although the free expansion assumption v — r/t is not applicable, v m now measures a true velocity 
of the photospheric radius (not only the matter flow speed, as in type IIP). 

• Now the original Baade's idea works for measuring the radius R p h by integrating di? p h = Vp^dt 
(of course, with due account of scattering, limb darkening/brightening etc. in a time-dependent 
modelling) . 




(11) 



H w 79Ii7 km/s/Mpc. 



(12) 
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• The observed flux then gives the distance D from the system ((8]). 

The constraining of cosmological parameters and our understanding of Dark Energy depend strongly on 
accurate measurements of distances in Universe. SNe Iln may be used for cosmology as primary distance 
indicators with the new DSM method. Application of EPM and SEAM requires crafting a best fitting 
hydro model for each individual SN. This procedure is in principle simpler in DSM. The case of SN 2006gy 
shows that the DSM distance agrees well with other most reliable techniques when the correct model is 
used, without the assumption on free expansion which is needed for EPM and SEAM. 
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